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ABSTRACT 

We study statistical intrinsic shape of star-forming BzK galaxies (sBzK galaxies) at z ~ 2 in both rest-frame 
UV and rest-frame optical wavelengths. The sBzK galaxies are selected down to Kab = 24.0 mag in the 
Great Observatories Origins Deep Survey-South (GOODS-S) and Subaru-XMM Deep Survey (SXDS) fields, 
where high-resolution images from Hubble Space Telescope are publicly available. 57% (583) of all 1028 
galaxies in GOODS-S show a single component in the Advanced Camera for Survey (ACS)/F850LP image. 
As Wide Field Camera (WFC3)/F160W images cover only some part of GOODS-S and SXDS, 724/1028 and 
2500/29835 sBzK galaxies in the GOODS-S and SXDS have the WFC3 coverage. 86% (626) and 82% (2044) 
of the sBzK galaxies in WFC3/F160W images appear as a single component in the GOODS-S and SXDS, 
respectively. Larger fraction of single-component objects in F850LP images represents multiple star-forming 
regions in galaxies, while they are not so obvious in the F160W image which appears smoother. Most of the 
single-component sBzK galaxies show Sersic indices of n ~ 0.5 — 2.5, in agreement with those of local disk 
galaxies. Their effective radii are 1.0 — 3.0 kpc and L5 — 4.0 kpc in F850LP and F160W images, respectively, 
regardless of the observed fields. Stellar surface mass density of the sBzK galaxies is also comparable to that 
of the local disk galaxies. However, the intrinsic shape of sBzK galaxies is not a round disk as seen in the 
local disk galaxies. By comparing apparent axial ratio (6/a) distributions of the sBzK galaxies with those 
by assuming tri-axial model with axes A > B > C, v^t found their intrinsic face-on B/A ratios peak at 
B /A = 0.70 and B/A = 0.77 — 0.79 in the rest-frame UV and optical, respectively and are statistically more 
bar-like than that of the local disk galaxies. The intrinsic edge-on C /A ratios in both rest-frame UV and optical 
wavelengths peak at 0.26, which is slightly larger than that of the local disk galaxies. Possible origins of this 
bar-like structure are bar instability, galaxy interaction, or continuous minor mergers. Study of galaxy structure 
evolution in cosmological simulations is desirable to examine the most likely scenario. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: structure 



1. INTRODUCTION 

Disk galaxy is one type of galaxies commonly seen in the 
present-day universe. Origin of the disk galaxy is, therefore, 
a crucial key to understand the galaxy formation and evolu- 
tion. Observational studies suggest that disk galaxies already 
exist at z ^ 1. Despite some luminosity evolution up to 
z ^ 1 (e.g.. IScarlata et aDl2007h . stella r surface mass den - 
sity does not show significant evolution (iBarden et al.ll2005h . 
Likewise, there is almost no significant evo lution of disk 
galax ies seen in the stellar mass function (e. g..lPannella et al] 
20091) and in the sc ale length function (e.g.. iLillv etalj|1998l: 
Sargent et al.l 120071) . These studies imply that the origin of 
disk galaxies is beyond z ^ \. However, obvious counter- 
part of disk galaxies has not been identified at z 3. Stud- 
ies of Lyman break galaxies (LBGs) at z 3 reveal that 
most of the LBGs show surface brightness distribution con- 
siste nt with Sersic profile with Sersic index n around one 
(e.g. iSteidel et al."1996t lAkivama etaI1l2008l) . Nevertheless, 
lAkivama et al. (2008y showed that the stellar surface mass 
density of the LBGs is too high to plausibly evolve into the lo- 
cal disk galaxies, suggesting that the LBGs at z ~ 3 are more 
likely to be progenitors of the local elliptical galaxies. In ad- 
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dition, large clustering amplitude of the bright LBGs indicates 
that they reside in massive dark matter halo and are seemingly 
not a direct progenitor of the present-day disk g alaxies (e.g., 
iGiavaUsco & Dickinsonll2001l:lOuchi et alj|200lh . 

At z ~ 2, galaxies can be efficiently selected by color se- 
lection (e.g.. ISteidel etalj r2004; Dad d i et al. 2004) . The clus- 
tering study of star-forming BzK (sBzK) galaxies, which ar e 
selected based on B ~ z and z — K color ( iDaddi et alj|2004t) . 
showed that faint sBzK galaxies {Kab < 23.2 mag) reside in 
dark matter halos with comparable mass to that of th e local 
disk g alaxies (|Havashi et al. 2007). Forster Schreibe r et alj 
(120091) studied kinematic properties of 80 star-forming galax- 
ies at 1.3 < z < 2.6 through spectroscopic imaging obser- 
vations of Ha emission and found that about 30% of their 
sample is rotationally supported with similar velocity-size re- 
lation to the present-day disk galaxies, albeit larger velocity 
dispersion. The results seem to imply that star-forming galax- 
ies at z 2 are likely to be disk galaxies with thicker disk 
than the local disk galaxies. Unfortunately, the intrinsic shape 
of these galaxies at z 2 had not been studied yet to confirm 
the presumption. 

The intrinsic shape of galaxies is related to distributions of 
their apparent axial ratios (6/a) and can be constrained by 
assuming a t ri -axial model with axes A > B > C (e.g., 
RvdenI pol iPadilla & StraussI l2008t lUnterborn & RvdenI 
20081) . I Yuma et alj (1201 IL hereafter paper I) studied the in- 



trinsic shape of sBzK galaxies in the GOODS-North field 
by using high-resolution imaging data obtained with the Ad- 
vanced Camera for Surveys (ACS) camera on Hubble Space 
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Telescope (HST). Among 1029 sBzK galaxies selected down 
to Kab — 24.0 mag, 54% of them showed a single com- 
ponent in ACS/F850LP images, of which the effective wave- 
length corresponds to the rest-frame UV 3000 A). Most of 
these single-component sBzK galaxies show surface bright- 
ness distribution with Sersic indices of n ~ 1, suggesting 
that they have disk-like structure. Besides, they found that 
the sBzK galaxies largely distribute in the same region of 
the size-mass diagram as z ^ — 1 disk galaxies, suggest- 
ing a comparable stellar surface mass density. However, pa- 
per I found that the sBzK galaxies at z ^ 2 apparently pre- 
fer a bar-like or oval shape rather than a round disk seen in 
the local disk galaxies. In case of the sBzK galaxies (paper 
I), the mean face-on axial ratio (B/A) is 0.61 and the mean 
disk thickness C/A is 0.28, while they ar e 0.90 - 0.92 and 
0.21 — 0.22, respectively, for the local disks jPadilla & StraussI 
120081: lUnterborn & Rvdenl2008l) . In other words, despite sim- 
ilar Sersic index and comparable stellar surface mass density, 
the sBzK galaxies at z 2 statistically show a bar-like or 
oval intrinsic shape different from disk galaxies in the local 
universe that show a round and flat disk shape. 

Although the conclusion about the intrinsic shape of sBzK 
galaxies is quite clear, some problems still remain for further 
investigations. The relevant problem is difference of effective 
wavelengths of the images used to derived the apparent b/a 
between galaxies at z ^ 2 and those at z 0. The apparent 
b/a of the sBzK galaxies in the GOODS-N field was deter- 
mined in the ACS/F850LP images corresponding to the rest- 
frame UV wavelength where the luminosity is mainly dom- 
inated by star-forming activity. In contrast, studies of local 
galaxies were carried out in the optical wavelength, which is 
dominated by more evolved stellar component. Thus their 
structures can possibly be different. Recent morphological 
studies of galaxies at z ^ 2 demonstrated that a galaxy may 
look different in different wavelengths (e.g., iCameron et al.l 
1201 U IConseUce et al.ll20TTl) . IConselice et al.1 (1201 ll) studied 
the stellar-mass selected galaxies at 1 < z < 3 and found that 
galaxies appearing peculiar in the rest-frame UV images often 
show compact morphology in the rest-frame optical images. 
Though this result is partially due to the resolution difference 
between images in two wavelength ranges, the morphologies 
of galaxies seem to be varied when observing in different 
wavelengths. Accordingly, it is important to study the intrin- 
sic shape of galaxies in the rest-frame optical wavelengths. In 
this paper, we examine the intrinsic shape of sBzK galaxies at 
z - 2 in GOODS-South and SXDS fields in both rest-frame 
UV and optical wavelengths. In addition to determine the dif- 
ference in structure between two different wavelength ranges, 
this work increases at least twice the number of galaxies used 
to derive the intrinsic shape in both wavelengths, which im- 
proves the statistical significance of the determined shapes. 

This paper is organized as follows. In Section|2] we provide 
an overview of the data sources and photometric catalogs of 
objects in both GOODS-South and SXDS fields. In Section 
[51 we describe how we determine the photometric redshifts 
and stellar properties of objects. Section |4] is given to a de- 
scription of our sample selection method and counterparts of 
sBzK galaxies in high-resolution images. Section |5] presents 
structural analysis and accuracy of structural parameters. In 
Section |6l we present results of structural parameters, stellar 
surface mass density, and intrinsic shape of the sBzK galaxies 
derived in both rest-frame U V and optical wavelengths as well 
as the comparison between different wavelengths. In Section 



|7] we compare our results with those of the local disk galaxies 
and discuss possible scenarios for the origin of the intrinsic 
shape of the sBzK galaxies. Finally, we summarize our re- 
sults in Section [8] We use a standard ACDM cosmology of 
Q„, = 0.3, ^ 0.7, and Hq = 70km s^^ Mpc^^ and AB 
magnitude system. 

2. DATA SOURCES AND PHOTOMETRIC CATALOGS 
2.1. GOODS-S 

Data in the GOODS-South field were collected as a part 
of the Great Observatories Origins Deep Survey (GOODSfl 
The optical HST/ACS data were taken from the v2.0 
data products of the G OODS HST/ACS treasury program 
dGiavalisco et al.ll2004l) . They consist of four band imaging 
data: F435W, F606W, F775W, andF850LP hereafter referred 
to as the Vgeo. *775- and zgso band, respectively. The 

images cover an area of ^ 160 arcmin^. The pixel scale 
of the images is 0."030 pixel^^. The 5cr limiting magni- 
tudes of the -8435, Vgeo, ^775, and zs5q images at l."0 di- 
ameter aperture are 27.9, 27.8, 27,3, and 27.0, respectively. 
The U and R data were obtained with yLT /VIMOS from the 
ESO/GOODS program (iNonino et alj|2009l) . The final U and 
R images in the area common to the GOODS HST/ACS ob- 
servations reach 5a limiting magnitudes of 28.61 and 28.06 
mag, respectively (at l."0 diameter aperture). 

The NIR VLT/ISAAC data in J, H, bands were taken 
from the GOODS/I SAAC data release final version (2.0) 
(lRetzlaffetal.1120101) . The pixel scale of the ISAAC images 
is 0."150, which is chosen so that one ISAAC pixel subtends 
exactly a block of 5 x 5 ACS pixels. The images in each band 
consist of at least 24 section^ over the GOODS-S field with 
different seeing size varying from 0."4 to 0."6. The 5cr limit- 
ing magnitudes for point sources of J = 25.0, H = 24.5, and 
Ks = 24.4 mag are reached within 75% of the survey area 
( lRetzlaffetal.1120101) . The mid-infrared (MIR) data were ob- 
tained from deep observations with the Infrared Array Cam- 
era (IRAC) on the Spitzer Space Telescope (SST). We used 
3.6fim and 4.5 fim images from both DRl and DR2 provided 
by the SST Legacy Science program. Limiting magnitudes 
of the 3.6/im and 4.5/im images are 24.3 and 24.4, respec- 
tively (5(7 at 2. "4 diameter aperture). The effective area of the 
GOODS-S field we study is - 160 arcmin^. 

In order to study the structural parameters of sBzK galax- 
ies in the rest-frame optical wavelength, the high-resolution 
NIR image is necessary. At z ^ 2, the rest-frame opti- 
cal wavelength 5000A) is redshifted to ~ 1.6/im, corre- 
sponding to the F160W band. We used the publicly available 
data from the Cosmic Assembly Ne ar-Infrared Deep Extra- 
galactic Legacy Su rvey (CANDELS: lKoekemoer et al.]|201 It 
iGrogin et al.ll201 ll) Q. This survey consists of five fields, one 
of which is the GOODS-S field. The data were taken from 
data release vO.5 of both deep JH and wide JH subregions, 
covering about 70% of the ACS/GOODS-S field with 10 ob- 
serving epochs. Limiting magnitude of the WFC3/F160W 
image (hereafter -ffigo) is ^ 25.6 — 26.0 depending on the 
observing epochs {5<j at l."0 diameter). The pixel scale is 
0."060 pixel- 1. 

2.2. SXDS 

* http ://www.stsci. edu/ science/goods/ 

' There ai'e two additional sections for J and Ks data. 

* http://candels.ucolick.org/index.html 
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The deep multi-wavelength data in the SXDS/UDS field 
are also publicly available. The optical Subaru/Supiime- 
Cam images {B, V, Rc, and z') were obtained from the 
Subaru-XMM Deep Survey (SXDS; iFurusawa et"an 120081) 
with the original pixel scale of 0."202. ba limiting magni- 
tudes at 2."0 diameter aperture are 27.7, 27.2, 27.1, 27.0, and 
26.0 mag for B, V, Rc, and z', respectively. The NIR 
data are fr om DR8 version of th e UKIDSS Ultra Deep Sur- 
vey (UDS; [Lawrence et al.ll2007l) . The J, H, and images 
with the original pixel scale of 0."268 were taken with WF- 
CAM on UKIRT. The limiting magnitudes are 24.9, 24.2, 
and 24.6 mag for J, H, and Kg, respectively (5(t at 2."0 diam- 
eter). The S'5r/IRAC images with the pixel scale of 0."600 
were obtained from the Spitzer Public Legacy Survey of the 
UKIDSS Ultra Deep Survey (SpUDS; PI: J. Dunlop). The 
limiting magnitudes are 24.8 and 24.5 for 3.6/im and 4.5/im, 
respectively. We used the overlapped area of these data, which 
effectively covers ^ 0.67 deg-^. The high-resolution NIR im- 
ages (iJi6o), which is used to study the galaxy structure, were 
taken from the CANDELS data release vl.O covering central 
~ 210 arcmin^ of the original SXDS. The Su limiting mag- 
nitude of the HiQQ image is 27.1 mag for point sources. It is 
noteworthy that the CANDELS images in the SXDS are ap- 
proximately one magnitude deeper than those we used in the 
GOODS-S field. 

2.3. Photometric Catalogs 

The photometric catalog was made by first selecting the 
objects in the ISAAC/Kg images and then making photom- 
etry in all images (except for WFC3 data) ranging from UV 
to MIR wavelength to construct spectral energy distribution 
(SED) of each object in the GOODS-S field. Due to seeing 
size difference among the ISAAC images, we first smoothed 
ISAAC/Ks images into one homogenous seeing size, i.e., 
0."6. The ACS/optical images were put into the same pixel 
scale as the ISAAC/NIR images by summing 5x5 pixels of 
the original images and preserving the fluxes. The images 
were then smoothed to have the FWHM common to the Kg- 
band images. Because of the dimension differences of images 
in different bands, objects were extracted separately above 
la w ith 5 connected pixels by SExtractor (iBertin & ArnoutsI 
119961) . For the ACS images, aperture photometry was carried 
out separately in each band centering at the position extracted 
in Kg band with 0."8 diameter aperture, which provides the 
best signal-to-noise (S/N) ratio for the PSF homogenized im- 
ages. The total magnitude of objects in ACS images was cal- 
culated from the aperture magnitude scaled up to the SExtrac- 
tor's MAG_AUTO in Kg band. For VIMOS and ISAAC images, 
the total magnitude was obtained directly from MAG_AUTO of 
each object. For IRAC images, the total magnitude was calcu- 
lated by using the aperture photometry at 2. "4 diameter and 
aperture correction. The correction factors were determined 
by generating the artificial objects with given total magnitudes 
and recovering their aperture photometry (paper I). 

The photometri c catalog of objects in SXDS is described by 
lYabe et al.1 jMl) . The NIR images are aligned and homog- 
enized to the optical images. Aperture photometry of objects 
in each image was made at the position extracted in the Kg- 
band image by using dual-image mode in SExtractor with an 
aperture of 2."0 diameter. Total magnitudes are calculated by 
applying the correction factor, which is determined by scal- 
ing the aperture magnitude to the MAG_AUTO in isTg-band im- 
age, to the aperture photometry in each band. For the IRAC 
images, aperture photometry was carried out at 2. "4 diameter 



aperture and was corrected to the total photometry in the same 
manner as done in the GOODS-S field. 

3. PHOTOMETRIC REDSHIFTS AND SED FITTING 

The photometric redshift (zphot) of each object in both 
GOODS- S and SXDS fields was determined by using 
Hyperz (iBolzonella et aljr2000h . The observed spectral en- 
ergy distributions (SEDs) of each individual object, which 
consist of photometry in UBVRizJHKg and IRAC 3.6 and 
4.5yum band, was fitted to the template SEDs covering E, 
SO, Sa-Sd, Im, and starburst types. The redshift range is 
set to be < Zphot < 6.0 with a step of Sz = 0.03. 
The obtained Zphot was compared with the a vailable spec- 
troscopic redshift catalogs in both GOODS-S ( iPopesso et alj 
120091: iBalestra et al.l 12010) and SXDS/UDS (Si mpson et al 
201 1, inprep.; Akiyama et al. 2011, in prep.; ISmail et al] 
l2008h . Although there are catastrophic redshift outliers, the 
Zphot is mostly consistent with the Zspec with uncertainties o-^ 
(standard deviation of Az/[1 + Zspcc]) of 0.05 and 0.03 in 
GOODS-S and SXDS/UDS, respectively. 

Stellar mass, star formation rate (SFR), age, and color 
excess {E(B-V)) of objects in the catalogs were obtained 
from SED fitting with the standard minimization. The 
SED fitting w as performed by using the SEDfit program 
(ISawickill20I2b . which is the evolved versio n of the SED fit- 
ting program used in (ISawicki & YeellI998h. The m a in dif- 
ferences are that the new version uses ICalzetti et alj ( 120001) 
law instead o f the ICa lzetti ( 199% and the population synthe- 
sis model by Bruzual & Charloj ( 120031) instead of those by 
iBruzual & Chariot (il993j) . BC03 population synthesis code 
(iBruzual & Charloll 120031) was used to construct the model 
SEDs with various star formation histories (instantaneous 
burst, exponential declining SF, constant SF ). The Salpeter 
initial mass function (IMF; ISalpeteiilI955l) with the mass 
range of 0.1 — 100 Mq and metallicity of I.QZq were as- 
sumedQ The dust attenuation law by ICalzetti et al.l (120001) 
was adopted with E{B - V) range of 0.0 - 1.0 mag (0.02 
mag step) and the redshift was fixed to the Zphot determined 
by Hyperz to reduce the number of free parameters. The un- 
certainty is obtained based on Monte Carlo simulation at 68% 
confidence level. The median uncertainties for derived stel- 
lar masses are Alogil/*/AfQ = 0.26 dex and 0.38 dex for 
samples in GOODS-S and SXDS, respectively. 

4. SAMPLE SELECTION 

4. 1 . Star-Forming BzK Galaxies 

Star-forming galaxies at z ^ 2 were selected by applying 
the BzK selection method ( iDaddi et al.l 2004) to the photo- 
metric catalogs described above down to Kg ~ 24.0 mag. 
Because the respons e functions used in t his paper are not iden- 
tical to those used by lDaddi et al J (120041) . the color corrections 
were determined in both GOO DS-S and SXDS field. The 
stellar spectra bv lPicklell (119981) were convolved with the re- 
sponse functions including the detector's quantum efficiency 
and atmospheric transmission. Both B — z and z — K colors 
derived from filter system used i n the GOODS-S ar e not sig- 
nificantly different from those by lDaddi et al.l (12004ft . Thus no 
correction is required for the GOODS-S field. In contrast, we 
adopted the color corrections for the SXDS field as follow: 

' Metallicity of most star-forming galaxies at 2 ~ 2 is slightly lower 
than solar abundance; ^ O.&Zp) or larger for ^ 10^^ Mp) ga laxies (e.g., 
IHavashi et aP [20091 : lYoshikawa etaH [20Tol : IWuvts et alj 12012) . Here we 
adopted the I.OZq model rather than the 0.2Zq model. 
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Figure 1. BzK diagram (-B435 — 2:^50 and zg^o — K) of sBzK galaxies in the GOODS-S field. The sBzK criterion is shown in a black dashed line. Left: Blue 
open triangles are the single-2:g5o sBzK galaxies. Right: Blue solid triangles represent those with multiple components in the original 2350 image. All objects 
with Ks < 24.0 mag are also shown with black dots. 




Figure 2. Normalized histograms of Ks-band magnitudes of sBzK galaxies in GOODS-S {left panel) and SXDS (right panel). The histograms are normalized 
by the number of all sBzK galaxies in the field. (For SXDS, the histograms are normalized by the number of all sBzK galaxies detected in the CANDELS 
field, which is part of the SXDS.) Left panel: Blue dashed and red dotted histograms represent the sBzK galaxies in the GOODS-S field appearing as a single 
component in the zgso and Hieo images, respectively. The black solid histogram is for the entire sBzK sample in the field. Right panel: The single-component 
sBzK galaxies in SXDS are illustrated with magenta dot-dashed histogram, whereas all sBzK galaxies appearing in the CANDELS field of SXDS are shown with 
black solid line. 



(B - ^;)Daddi+04 = (-B - z)sxDS + 0.3, (1) 
(.Z - i4:)Daddi+04 = - .f'^)sXDS + 0.1. (2) 

Then the original sBzK criterion, BzK = (z — A') — {B — 
z) > —0.2, was applied to construct the sBzK samples in both 
fields. The candidates with Zgpec < 1-4 and Zgpec > 2.5 were 
excluded from the samples. Eventually, we have 1028 and 
29835 sBzK galaxies in GOODS-S and SXDS, respectively. 

4.2. Cross-Identification with HST Images 
4.2.1. ACS/F850LP 

In order to verify the intrinsic shape of sBzK galaxies in 
the same rest-frame UV wavelength as in the GOODS-N field 
(paper I), we used the zs5o images that are publicly avail- 
able in the GOODS-S field. In section |2l we constructed 
the photometric catalog of objects selected in A's-band im- 
ages by using the smoothed ACS images. However, we also 
need to identify galaxies using high-resolution images in or- 
der to study structure of galaxies. The sBzK galaxies are re- 
identified in the original zgso image, which corresponds to 
rest-frame ~ 3000A. An object in the zgso images is identi- 
fied as a BzK counterpart if it is detected above 2a threshold 



with 5 connected pixels in SExtractor and locates within l."0 
radius centered at the sBzK position in Kg band. Similar to 
paper I, the radius was chosen as it is twice of the median 
1/2FWHM of the samples in Kg images. The sBzK galax- 
ies are divided into 2 groups; the single-component and the 
multiple-component sBzKs, depending on the number of their 
counterparts within the radius. 583 (57%) sources are identi- 
fied as single-component sBzK galaxies, while the multiple- 
component sample contains 438 galaxies including six objects 
(3 pairs) of (interacting) sBzK galaxies. Note that the ratio 
of single-component sBzK galaxies is in agreement with that 
(54%) in GOODS-N (paper I). The contamination rate of for- 
ground/background galaxies within l."0 radius is estimated 
to be 14% for the zs5q image. As we select the single- and 
multiple-component galaxies in both zsso and iJieo images, 
we explicitly denote the single-component sBzK galaxies in 
the Z850 images as the single-zgso galaxies and the multiple- 
component sBzK galaxies in the 2359 images as the multiple- 
Z8.50 galaxies to avoid any confusion. Figure [T] shows BzK 
diagrams of the single-zgso and multiple-zgso sBzK galaxies 
in the z^^o images. Distributions of both single- and multiple- 
zgso sBzK galaxies are almost the same. They show simi- 
lar distribution to the whole sBzK sample in the GOODS-S 
field as well. A histogram of Kg magnitude of the single- 
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Single-comp sBzK (GOODS-S),. 




Multiple-comp sBzK (GOODS-S)' 



2 4 
B-z 



2 4 
B-z 



Single-comp sBzK (SXDS) 




Multiple-comp sBzK (SXDS) 




Figure 3. BzK diagram of the sBzK galaxies detected in the f/ie o image in GO ODS-S (top panel) and SXDS {bottom panel). The B — z and z — K colors for 
objects in SXDS were corrected to match the filter system used by lDaddi et alj {20041) (see text for more details). The sBzK criterion is shown in a black dashed 
fine, top panel: Red open triangles are the single-i/ieo sBzK galaxies in the GOODS-S field, while red solid triangles represent the multiple- Hieo galaxies. 
bottom panel: The single- and multiple- //leo sBzK galaxies in SXDS are shown with magenta open and solid circles, respectively. All objects with Ks < 24.0 
mag are also shown with black dots. The objects that are above the sBzK criterion but are neither classified as single-component nor multiple-component objects 
are those outside the Higo images, which cover a part of the entire fields of GOODS-S and SXDS. 



zgso sBzK galaxies is shown in the left panel of Figure |2] 
together with the histograms of all sBzK galaxies. The figure 
shows that the single-Z85o galaxies similarly distribute in the 
magnitude range of 21 — 24 mag as compared to the total sam- 
ple in the field. The total number of the sBzK galaxies that are 
able to be identified to either single or multiple components 
in the high-resolution zs5q images is 1021 out of all 1028 ob- 
jects in the GOODS-S field. The remaining seven candidates 
are too faint to be detected in the original zq^q images, though 
they were detected in the iiTs-band homogenized zs5o images, 
which were used in the sBzK selection process. 

4.2.2. WFC3/F160W 

The high-resolution WFC3/F160W (-ffieo) images are used 
in order to determine the structure of the sBzK galaxies in 
the rest-frame optical wavelength range. The -ffieo images 
have the effective wavelength corresponding to the rest-frame 
~ 5300A at z ^ 2. We applied the definition identical 
to that used for the ACS/F850LP images in order to clas- 
sify the single or multiple components, i.e., the sBzK coun- 
terparts in the high-resolution images are identified if be- 
ing detected above 2<t and 5 connected pixels within l."0 
radius from the sBzK center. 626 (86%) single-component 
and 98 multiple-component sBzK galaxies are obtained in the 
GOODS-S, and we have 2044 (82%) single-component and 
456 multiple-component galaxies in the SXDS. Note that be- 
cause the WFC3/F160W data are not available for the whole 
area of the GOODS-S and SXDS fields, only some fractions 



of the sBzK galaxies are used in this paper. The contamination 
rate for the -ffieo images is 7% within l."0 radius. Although 
most of sBzK galaxies located in the -ffieo images (> 80%) 
are classified as a single component in both fields, the single- 
component fraction somehow does not agree to each other 
within Icr uncertaint>0, i.e., 86 ± 1% for GOODS-S and 82 ± 
1% for SXDS. However, they are still in agreement with each 
other within 2<t range. We refer to the single- and multiple- 
component sBzK galaxies in the -ffieo images as the single- 
HiQo and the multiple-iJieo sBzK galaxies, respectively. The 
BzK diagrams of both single- and multiple-ifieo sBzK galax- 
ies in both fields are illustrated together with the whole sBzK 
sample in Figure 13 The single-iJieo sBzK galaxies in both 
fields seem to show similar distributions in the BzK diagram 
as compared to the entire sample of sBzK galaxies. A'^-band 
magnitude histograms of the single-TJigo sBzK galaxies in 
GOODS-S and SXDS are shown in the left panel (red dot- 
ted) and right panel (magenta dot-dashed) of the Figure |2] 
respectively. Histograms of the single-iJieo galaxies in both 
fields are similar to those of the whole sBzK galaxies shown in 
the black solid lines. These imply that the single-iJigo sBzK 
galaxies can be good representatives for the entire populations 
at least in term of their colors and rest-frame optical magni- 
tudes. 



cr^ = (1 — /) * f /N, where / denotes the fraction and A'^ is for number 
of the sample. 
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Table 1 

Number of sBzK galaxies detected in high-resolution images 







GOODS-S 




SXDS 


Type of sBzK 


^850" 


Hl6o-smoothed zgso 




^^160 


Single-component 
Multiple-component 


583 (57 ± 2%) 
438 (43 ± 2%) 


718 (70 ±2%) 
310 (30 ±3%) 


626 (86 ± 1%) 
98 (14 ± 4%) 


2044 (82 ± 1%) 
456 (18 ± 2%) 



^ The number of all sBzK galaxies detected in the original zgrM images is 102 1 , which is not identical to 
the number of all sBzK in GOODS-S field (i.e., 1028). This is because the sBzK galaxies were selected 
based on the zgso images homogenized to match PSF of the Ks-hwd images instead of the original 
images (see text). 



4.2.3. High Fraction of Single-Component sBzK Galaxies in 
WFC3/F160W 

As described above, the sBzK galaxies are categorized as 
single-component or multiple-component galaxies according 
to their appearance in the high-resolution images both in the 
rest-frame UV (zgso) and the rest-frame optical (i/ieo) wave- 
lengths. In the GOODS-S field where high-resolution images 
are available in both rest-frame UV and optical wavelength, 
57% of all sBzK galaxies are identified as single-Z8.5o galax- 
ies, while 86% of the galaxies located in the i/igo images are 
the single-i/igo samples. There are two plausible reasons re- 
sponsible for the difference in fraction of single-component 
objects in different wavelengths. The first reason is the dif- 
ference in the resolution of images between in zgso and in 
iJieo bands. The zs5o images have about two times better res- 
olution than the i^igo images; FWHM of the PSF in z^rM and 
HiQQ images are ^0."1 and ~0."18, respectively. The second 
reason is that the difference is real. The rest-frame UV com- 
ponent of a galaxy is mainly dominated by star-forming ac- 
tivity, whereas the relatively low-mass, long-age stellar com- 
ponent tends to dominate the rest-frame optical wavelengths. 
Consequently, components of one galaxy seen in the rest- 
frame UV (i.e., zgso) are possibly just the star-forming regions 
in the galaxy, not the entire galaxy. 

In order to examine which reason is responsible for this 
different fraction, we first smoothed the original zgso image 
to the HiQo resolution. Then we re-detected and re-selected 
the single-component and multiple-component sBzK galaxies 
from the i^ieo-smoothed Z85o images. The number and per- 
centage of the single- and multiple-component sBzK galaxies 
are summarized in Table [T] As seen in the table, fraction of 
the single-zgso galaxies increases significantly after homog- 
enizing the Z8,5o images to match the PSF of the Hiqq im- 
ages; however, the fraction is still less than that for the single- 
Hieo sample. The results indicate that the difference in ap- 
pearance of sBzK galaxies between in the rest-frame UV and 
optical wavelength is partially due to the difference in resolu- 
tion of images obtained with ACS and WFC3. On the other 
hand, more than half of the single-i/igo galaxies seen as mul- 
tiple components in the Z85o images indicate real difference 
between the morphology of galaxies in the rest-frame UV and 
that in the rest-frame optical wavelength. 

5. STRUCTURAL ANALYSIS 

We carried out the structural analysis in the same process as 
done in paper I. Briefly, the structural parameters of the sBzK 
galaxies are obtained by fitting their two-dimensional surface 
brightness dist ributions in zsbo or ffiso images with a single 
Sersic profile ('Sersic" 1 9631: ISersidl 1 9681) using GALFIT ver- 
sion 3.0 ( iPeng et al...2010l) . The Sersic profile is formulated 



Table 2 

Number of the single-component sBzK galaxies before and after 
GALFIT analysis 





Number of sBzK galaxies (After/before GALFIT)" 


Observed Fields 


Single-2850 Single-Z/ieo 


GOODS-S 


467/573 (82%)'' 550/626 (88%) 


SXDS 


1878/2044 (92%) 



" The objects before GALFIT are those used in GALFIT process. The 
objects after GALFIT are those with < 1-0 and errors within 30% of 
the derived parameter values. 

Although there are 583 single-^gso sBzK galaxies, 10 of them are at the 
image edge which cannot be used to study the structural parameters. 
" The ACS/F850LP images are not available in the SXDS field. 

by I{r) = /e cxp { — fc„[(^)^/" — 1]}, where /g is surface 
brightness at an effective radius r,; and fc„ is a variable pa- 
rameter coupled with Sersic index n so that covers half of 
the total flux. The de Vaucouleurs and exponential profiles are 
special cases of the Sersic profile when n = 4 and n = 1, re- 
spectively. Seven free parameters of each object, i.e., position 
(x, y), magnitude, n, r^, axial ratio (b/a), and position angle 
(PA), are determined by minimization. Initial values of all 
parameters are given from SExtractor outputs, except for the 
Sersic index which is initially set at n = 1.5. As explained 
in paper I, the GALFIT results do not significantly depend on 
the initial Sersic index. Sigma images required for calculating 
the and fitting errors are constructed from variance maps 
produced during the reduction process. 

The analysis was executed for all single-component sBzK 
samples, i.e., the single-Z85o & single-iJigo sBzK galaxies 
in the GOODS-S and the single-i/igo sBzK galaxies in the 
SXDS. We excluded objects with xt larger than one and con- 
taining inconvincibly large fitting uncertainties, i.e., objects 
with parameter errors more than 30% of the parameter val- 
ues. The excluded galaxies are mostly too faint to obtain re- 
liable results. The resultant numbers of sBzK galaxies left 
after GALFIT analysis are listed in Table |2l Over 80% of 
the single-component sBzK galaxies in both fields show good 
fitting results of GALFIT. Examples of the GALFIT results 
are demonstrated in Figures HE] and m for the 

single-2850j 

single-i/i6osBzK galaxies in the GOODS-S and the single- 
HiQo sBzK in the SXDS, respectively. The figures show 
the observed images, best-fitting models, and residuals, along 
with one-dimensional light profile visualizing the agreement 
between the observed and modeled profiles. 

5.1. Accuracy in Determining Structural Parameters 
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z_spec=1.4671; 23.34 mag; 
Re=2.51kpc: n=0. 74; b/a=0.91 




z_spec=1 .8993; 24.32 mag; 
Re=1.77kpc; n=0.96; b/a=0.58 



z_spec=2.2865; 24.69 mag; 
Re=1.02kpc; n=2.21; b/a=0.56 



z_spec=1 .8778; 24.72 mag; 
Re=2.23kpc; n=0.98; b/a=0.63 



z_spec=2.4799; 25.01 mag; 
Re=2.13kpc; n=1.42; b/a=0.55 
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Figure 4. Examples of 2D surface brightness modeling with GALFIT for the single-zgso sBzK galaxies in the GOODS-S field. The observed zg^o image, the 
best-fitting model constructed by GALFIT, and the residual (the observed image subtracted by the model) are aiTanged from the leftmost panel. The residual 
images are shown with the narrower scale range to clarify the residual part. North is up; east is to the left. The size of each panel is 2."0x2."0 coiTesponding 
to ~ 17 X 17 kpc at z ~ 2. The rightmost panel shows the azimuthally averaged ID surface brightness profile (black open circles with errors) with the profile 
of best-fitting inodel (red solid line) and normalized PSF (blue dashed line). The surface biightness of the residual is also shown in the bottom of the rightmost 
panel. 



In paper I, we carried out the Monte Carlo simulations to 
examine the accuracy of the derived structural parameters in 
the zgso images in GOODS-N field and found that the estima- 
tion of Sersic index becomes less reliable at Z850 > 24.5 mag. 
The effective radii derived for disk profiles are robustly recov- 
ered down to Z850 = 26.0 mag and more accurate than those 
in the spheroid case. Since the 2^50 images in the GOODS-S 
are 0.5 mag deeper than those in the GOODS-N, we conser- 



vatively assume the same accuracy as estimated in paper I. 

In case of the i^ieo images, we separately generated ^ 
8000 artificial objects in both GOODS-S and SXDS with ran- 
dom magnitude, effective radius, ellipticities, and position an- 
gles of 20 - 25 mag, 0."1-1."0, 0.1 - 1.0, and - 180°, 
respectively. The ranges of the input parameters were chosen 
to be the same as those of the sBzK galaxies. The artificial 
objects, which are equally constructed for disks (71 = 1) and 
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Figure 5. Same as Figure|5]but for the Hiso images of the single-Zfieo sBzK galaxies in the GOODS-S field. Note that the slightly differences among PSFs 
shown here are due to the differences of adopted stars in the image taken at each epoch. 



spheroids (n ~ 4), were convolved with PSF of Hiqo im- 
ages and inserted into the original iJieo images. Then we 
performed the GALFIT and selected the good fits in the iden- 
tical way to that used for the sBzK galaxies; i.e., xt < 1-0 
and errors < 30%. The percentages of successful fittings 
for the Hiea artificial objects in both GOODS-S (86%) and 
SXDS (92%) fields are very similar to those for the single- 
i?i6o sBzK galaxies (Table|2]i. 

Figure |7] shows the accuracy in determining the structural 
parameters from Monte Carlo simulations in the GOODS- 
S field by GALFIT program along with Icr distributions of 
the recovered structural parameters. The left panel shows the 



Sersic index accuracy of both disk-like (n = 1; green solid 
square) and spheroid-like {n = 4; orange open square) ar- 
tificial objects. It is shown that the recovered Sersic indices 
are systematically underestimated in most cases. It is not sig- 
nificant in case of simulated disks; the median difference be- 
tween the input and recovered Sersic indices is no more than 
10% at all magnitude ranges. In contrast, the Sersic index es- 
timation for the simulated spheroids becomes less accurate at 
fainter magnitudes; their recovered Sersic indices go below 
n < 2.5 at i^ieo > 24.0 mag. In other words, we are able 
to distinguish the disk-like sample from the spheroid based 
on the Sersic index criterion n — 2.5 down to i/ieo < 24 
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Figure 6. Same as Figure|5]but for the -ff igo images of the single-Higo sBzK galaxies in the SXDS field. 
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mag in the GOODS-S. Similar to the simulations made in the 
zg5o images (paper I), the effective radius of a disk profile is 
better recovered than that of the spheroid one. The middle 
panel shows that the effective radii of the simulated disks are 
well recovered down to the faint magnitudes; however, they 
are systematically underestimated in case of spheroid about 
10 — 20% except for the faintest magnitude bin. However, 
this underestimation is not severe as we show in the next sec- 
tion that most of our samples are disk-like, of which the ef- 
fective radii can be accurately recovered. The other important 
parameter that we obtain from GALFIT is the apparent axial 
ratio (b/a), which will be used later in determining the intrin- 



sic shape of the galaxies. The resulting accuracy in determin- 
ing the b/a is shown in the right panel of Figure [T] The b/a 
ratios are robustly recovered in both types of profiles down 
to HiQQ = 24 mag, at which the recovered b/a begins to be 
underestimated. 

Similarly, the accuracy in recovering the structural param- 
eters of the artificial objects inserted in the SXDS is sum- 
marized in Figure [8] The simulations performed in both 
GOODS-S and SXDS fields show the similar accuracy in 
all parameters over most of the magnitude range. However, 
as the images in the SXDS are about one magnitude deeper 
than those in the GOODS-S, the accuracy calculated at the 
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Figure 7. Accuracy in determining structural parameters by GALFIT based on Monte Carlo simulations plotted as a function of input magnitudes for the 
GOODS-S field. From left to right panels, the figure shows the accuracy of Sersic index (n), effective radius (re), and axial ratio (b/a). Green solid squares 
represent the accuracy calculated for the artificial objects with an exponential profile (n = 1), while orange open squares show the accuracy for those with a de 
Vaucouleurs profile (n = 4). Error bars correspond to Icr of distributions. The horizontal dashed line at (riout — nin)/ni„ = —0.375 in the left panel indicates 
where a spheroid can be misclassified by GALFIT as a disk. 
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Figure 8. Same as Figure|7]but for the SXDS field. 
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faintest magnitude bins (i.e., Hiqo > 24.0 mag) in the SXDS 
is roughly better than that in the GOODS-S. The Sersic in- 
dices (left panel) are generally underestimated in both disk 
and spheroid profiles, but negligible (< 10.0%) in the disk 
case. The figure shows that the disk-like and spheroid pro- 
files are accurately distinguishable down to i^igo < 24.0 mag. 
Similar to the GOODS-S, the middle panel indicates that the 
effective radii are robustly derived in the disk profile, while 
they are averagely 10 — 20% underestimated in the case of 
spheroids. The axial ratios b/a are accurately determined in 
both profiles over the magnitude range as seen in the right 
panel of Figure[8] 

Surface brightness of galaxies can be different due to their 
inclination; edge-on galaxies tend to show higher surface 
brightness than the face-on galaxies at the same total mag- 
nitude. The sample selection can possibly be biased toward 
edge-on galaxies, which consequently results in lower a peak 
of b/a distribution. This selection bias was already tested in 
paper 1 by generating 10, 000 artificial objects with the 
random inclination angles (axial ratios) and uniformly dis- 
tributed magnitudes of 21 — 24 mag in the A', -band image. 
Their effective radii were assumed to be in the same range as 
the observed galaxies selected from the image at each magni- 
tude bin. They found that the detection rate remains constant 
above 95% regardless of inclination angle or axial ratio down 
to Kg = 23.0 mag (Figure 10 in paper 1). For the faintest mag- 
nitude bin (23 — 24 mag), the detection rate reaches 100% at 



b/a = 0.1 — 0.2 (edge-on case) and then becomes constant at 
^ 90% from b/a = 0.3 to 6/a = 1.0, which is a face-on case. 
Therefore, there seems to be no significant bias toward select- 
ing mostly edge-on galaxies. As we use the same Kg < 24 
mag limit as that in paper 1, we assume that brighter surface 
brightness of edge-on galaxies does not cause any bias on our 
sample selection. 

6. RESULTS 

6.1. Structural Parameters 

6.1.1. GOODS-S 

In this section, we present results of GALFIT of the single- 
Z850 and the single-iJigo sBzK galaxies in the GOODS-S 
field. Distributions of the derived Sersic indices and effec- 
tive radii are shown in Figure |9l The effective radii are in- 
dicated in the unit of kpc calculated by assuming the photo- 
metric redshift or spectroscopic redshift (if available) of each 
object. In the figure, we also show the histograms for the 
brighter sample; galaxies with zgso < 24.5 mag (top pan- 
els) and Hieo < 24.0 mag (bottom panels). It is seen in 
Figure |9] (a) that most of the single-zgso sBzK galaxies have 
n ^ 1 regardless of the magnitude limit. The median n for 
all and the brighter single-zgso sBzK galaxies are the same at 
71 ~ 1.17. Figure |9](b) shows that the effective radii of most 
single-Z85o galaxies are in the range of = 1.0 — 3.0 kpc, 
regardless of magnitude limit. The median are 1.89 and 
1.91 kpc for all single-^;85o galaxies and the brighter samples, 
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Figure 9. Histograms of the derived Sersic indices (n) and effective radii (re) of tlie single-^gso (top panels) and the single-Hiso (bottom panels) sBzK 
galaxies in the GOODS-S tield. The soHd histograms represent all single-^gso (blue) and single- /fiso (red) sBzK galaxies, whereas the open histograms are for 
the galaxies with zgso < 24.5 mag and Hiao < 24.0 mag for the single-zgso and the single- Hieo sBzK galaxies, respectively. The effective radii are described 
in unit of kpc calculated by assuming the photometric redshift or spectroscopic redshift (if available) for each object. 



respectively. 

Figures |9] (c) and (d) show the n and histograms for the 
single-i7i6o sBzK galaxies. Similar to those in the ^§50 im- 
ages, the Sersic indices of the single-iJieo sBzK galaxies 
mostly distribute around 71 ~ 1 with the median values of 
n ~ 1.09 and n = 1.08 for all samples and those with 
iJieo < 24.0 mag, respectively (Figure |9](c)). Figure |9](d) 
shows that the effective radii of the single-i7i6o galaxies are 
typically around = 1.5 — 4.0 kpc. The median for all 
single-ifigo sBzK galaxies is 2.52 kpc, similar to the median 
of fe = 2.64 kpc for the brighter sample. 

We performed Kolmogorov-Smirnov (K-S) tests to the n 
and distributions of all and the brighter sample and found 
that we could not reject the null hypothesis that the n and 
histograms of both samples are drawn from the same popu- 
lations at more than 98% confidence levels. We need to use 
as large number as possible of the sBzK galaxies in order to 
determine their intrinsic shape with sufficient statistical accu- 
racy. Since the distributions of both n and are consistent 
regardless of the magnitude limit, we keep using the whole 
samples of both single-zgso and single-i/igo sBzK galaxies 
to determine the intrinsic shape in section |7] We divide the 
sBzK galaxies into disk-like and spheroid-like samples ac- 
cording to their Sersic index. The disk-like galaxies are those 
with Sersic indices of 0.5 < n < 2.5, whereas the spheroid 
ones are those with n > 2.5. In both cases of the single- 
Z850 and single-iJieo sBzK galaxies, majority of them (72%) 
are similarly classified as disk-like galaxies. 

6.1.2. SXDS 



We illustrate the histograms of the Sersic indices (n) and 
the effective radii (rg) of the single-iJieo sBzK galaxies in 
Figures [To] (a) and (b), respectively. Figure [TO] (a) shows that 
the single-i/i6o sBzK galaxies in the SXDS mainly have the 
Sersic indices of n ~ 1 in both all and the brighter sample. 
The median values of both samples (n = 1.23) are very close 
to those obtained in GOODS-S. Using the same criteria as 
those in the GOODS-S, we categorized 67% of our all single- 
iJieo sBzK galaxies in the SXDS as the disk-like sample. Fig- 
ure [To] (b) also shows that the effective radii range between 
Te = 1.5 and = 4.0 kpc in both cases. The median r^. 
of the whole sample is 2.62 kpc, which is similar to that of 
iJieo < 24.0 sBzK galaxies (median = 2.76 kpc) and al- 
most equal to those in the GOODS-S (median rg = 2.52 kpc). 

6.1.3. Comparison between the Rest-frame UV and Optical 
Wavelengths 

As mentioned in section [T] a star-forming galaxy may have 
different structure when being observed in different wave- 
lengths. Star-formation activity probably dominates the rest- 
frame UV structure; therefore, studying the rest-frame optical 
structure is crucial to understand the galaxy morphology in 
the wavelength where more evolved stellar component dom- 
inates. In this section, we compare the structural parameters 
obtain from GALFIT between the rest-frame UV and optical 
wavelengths. Comparison of Sersic indices and effective radii 
of the single-component sBzK galaxies between in rest-frame 
UV and optical wavelengths in the GOODS-S field is shown 
in Figure [TT] Top left panel shows the Sersic index compar- 
ison for the galaxies detected in both zs5o and i/ieo images. 
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Figure 10. Histograms of the derived Sersic index (n) and effective radius (re) of tlie single-Hiso sBzK galaxies in the SXDS field. Solid magenta and open 
violet histograms show samples of all and -ffieo < 24.0 mag single-_ffi6o sBzK galaxies in the field, respectively. 
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Figure 11. Top panels: Comparison of Sersic indices (left) and effective radii (right) of the single-component sBzK galaxies in GOODS-S detected in both 
Z850 and Hie,a images. Bottom panels: Histograms of Sersic indices (left) and effective radii (right) of all single-component sBzK galaxies in GOODS-S. Blue 
and red histograms indicate the Sersic indices or the effective radii of the galaxies derived from 2:350 and i/igo images, respectively. Distributions are normalized 
to peak at unity. The median values are also shown by arrows with corresponding colors. 



It is important to note that there are 296 galaxies that are de- 
tected as a single component in both images. Sersic indices 
derived from both wavelength ranges of most galaxies are in 
agreement with each other la scatter of Sersic index ratio 
(nopt/nuv) ±0.2. Similar comparison but for the effective 
radius is shown in the top right panel. The effective radii de- 
rived in both wavelength range seem to be almost equal with 
little scatter. 

Bottom panels show histograms of all single-component 
sBzK galaxies detected in either zgso or iJieo images. The 
bottom panel shows that the Sersic index distributions of the 



sBzK galaxies in both wavelengths are averagely very similar 
and their median values are comparable with slightly larger 
Sersic index in the rest-frame UV wavelength. The median 
Sersic indices are 1.17 and 1 .09 in the rest-frame UV and rest- 
frame optical wavelengths, respectively. Histogram of the ef- 
fective radii is shown in the bottom right panel. The effective 
radii of all single-zg5o galaxies distribute in narrower range 
and are apparently smaller than those of all single-iJigo sam- 
ples. Their median ?> values are 1.89 and 2.52 in the rest- 
frame UV and optical wavelengths, respectively. The ratio 
of median in rest-frame optical to that in rest-frame UV is 
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Grids of Intrinsic Model Parameters 



9 10 
logCM/Mp) 

Figure 12. Size-mass relation of tlie single-i^iso sBzK galaxies with 0.5 < 
n < 2.5 both in GOODS-S (red open circles) and in SXDS (black open 
squares). The average uncertainties of the stellar mass and effective radius are 
shown in the top left corner of the figure with con'esponding colors. Dashed 
and solid contours illustrate size-mass distributions of disk galaxies at z = 
and 2 = 1, respectively (Harden etal. 2003). The black solid line is an 
average suiface mas s density for disk galaxies at 2 = 0—1. (log Ejv/ = 8.50 
with q = 0.5; iBarden et al. 2(305). The green line show s the typical size - 
mass relation of the elliptical galaxies in the local universe dShen et al.l2003n . 
Note that the stellar masses shown in the figure are all connected to Salpeter 
IMF. 

1.33, which agrees well with the ratio of 1 .37 found in the star- 
for ming galaxies at the sim ilar redshifts (BM/BX and LBGs) 
by ISwinbank et al.l (l2010h . The effective radii in the rest- 
frame optical wavelength are also roughly in agreement with 
those deri ved for the sBzK and BM/ BX galaxies at the s ame 
redshift bv lOverzier et al.l (120101) and lMosleh et al.l (120111) . re- 
spectively. It suggests the similar structures between the sBzK 
galaxies and the BM/BX or LBGs at z 2. 

The ratio of median seems to contradict with the re com- 
parison for individual galaxy that appears as a single compo- 
nent in both wavelength range (57% of all single-ifieo sBzK; 
top right panel) with the median re,opt/re,uv of 1.0. This is 
mainly because the galaxies with larger effective radii in the 
HiQo image appear as multiple components in the zs5o im- 
age; ^ 38% of all single-i^ieo sBzK galaxies are those clas- 
sified as multiple-Z85o galaxies and almost all of them have 
Te.opt larger than the median Ve^uv- It is worth to recall 
that higher fraction of single-iJigo sBzK against the single- 
Z850 sBzK galaxies are partly due to the poorer resolution in 
the WFC3/F160W image as compared to the ACS/F850LP. 

6.2. Stellar Surface Mass Density 

Figure [T2] shows the size-mass relation of the single- 
Hieo sBzK galaxies with 0.5 < n < 2.5 in both GOODS- 
S and SXDS fields and local galaxies. Because we derived 
the effective radii of the sBzK galaxies in the i/ieo images 
(i.e., the rest-frame optical images), the size-mass relation of 
the sBzK galaxies is able to be directly compared with the 
relation of galaxies in the local universe, where the effec- 
tive radii were measured in optical wavelengths. The stel- 
lar masses of local galaxies were corrected into Salpeter IMF 
with log{M Salpeter) - l oRiMKmupn) ^ 0-15 - 0.2 dex de- 
pendi ng onU — V colors ( lAkivama et al.l2008tlRudnick et alj 
l2003h . As seen in the figure, most of the sBzK galaxies with 
0.5 < n < 2.5 locate in the same region as the local disk and 
z ^ 1 disk galaxies locate. The result is similar to what we 
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found for the single-component sBzK galaxies in GOODS-N 
(paper I), confirming that the sBzK galaxies at z ~ 2 have 
comparable stellar surface mass density to the disk galaxies at 
z — 1. However, there are some outskirts that show larger 
surface mass density, which are closer to the typical value s 
of the local elliptical galaxies (green line; IShen et aLllIOOl . 

10*^ M„ are 



Note that the galaxies with stellar mass of 10^ 
those with the estimated photoz lower than 1 .4. 



6.3. Axial Ratio Distribution and Intrinsic Shape of sBzK 
Galaxies 

Although we found that the sBzK galaxies in both wave- 
lengths show the Sersic indices around one and have stellar 
surface mass density comparable to disk galaxies at z ~ — 1, 
it is still important to see whether or not their structure is a 
round disk as seen in the local disk galaxies. Since the in- 
trinsic shape of galaxies can be statistically examined from 
their observed distributions of apparent axial ratios (b/a), we 
show the b/a distribution and the resulting intrinsic shape of 
the sBzK galaxies in this section. 

6.3.1. Rest-Frame Win GOODS-S 

Histograms of the apparent b/a of the single-zgso galaxies 
with 0.5 < n < 2.5 in GOODS-S and GOODS-N fields (pa- 
per I) are shown in the left panel of figure [jj] Distribution of 
local disk galaxies by (iPadilla & Straussll2008L ~ 300, 000 in 
SDSS) is also shown. The distribution of samples in GOODS- 
S peaks at b/a ^ 0.35 and declines significantly toward 
b/a — 1.0. This trend agrees with that found in GOODS- 
N. We performed the K-S test and found that b/a distribution 
of the single-Z85o sBzK galaxies are drawn from the same dis- 
tribution with 92% confidence level. However, they are obvi- 
ously different from the distribution of the local disk, which 
is relatively flat in the range of 6/a 0.2 — 0.8. It is im- 
plied that structure of the sBzK galaxies in the rest-frame UV 
differs from that of the local disks. 

In order to study the intrinsic shape of galaxies, a tri-axial 
model is adopted with axes A > B > C. The face-on ellip- 
ticity (e = 1 — B/A) is assumed to be described by lognormal 
distribution with mean fi and dispersion a, while the edge-on 
thickness (C/A) is described by Gaussian distribution with 
mean pi^ and standard deviation cr^. One set of these 4 pa- 
rameters provides one apparent-6/a distribution after repeat- 
edly computing the apparent b/a at various random viewing 
angles. The parameters /i, a, /i-y, and were constrained 
by comparing the model distributions of 6/a with those of 
the observed samples. The ranges and steps of the param- 
eters are summarized in Table [3] The best-fitting parameter 
set was obtained by minimization method. The observed 
b/a distribution was first compared with the above parameter 
sets to obtain the primary result. Then we refitted the ob- 
served distribution with parameter set with finer grid (at least 
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Figure 13. Left panel: Normalized distributions of the apparent axial ratio (6/a) of the single-^gso sBzK galaxies in the GOODS-N (black dashed) and GOODS- 
S (blu e dot-dashed) fields. The best-fitting models are shown with solid lines with con'esponding colors. Distribution of local disk galaxies by Padilla & Straus^ 
120081) is also shown with thin solid histogram. Ri^ht panel: Distributions of the best-fitting intrinsic axial ratio: BjA distribution on the top right panel and 
C I A distribution on the bottom right panel with colors identical to the left panel. Y-axis is in an arbitrary unit. 



half of the current step) centered at the primary best-fitting 
result. The fitting uncertainties are calculated at 68% con- 
fidence level based on Monte Carlo realization. We derived 
the best-fitting set of parameters and repeated the process 500 
times by varying the observed distribution of b/a ratios within 
their Poisson noises. 

The final best-fitting results of the single-zgso sBzK galax- 
ies in GOODS-S are summarized in Table|4] For comparison 
purpose, we also show the results of the single-zggo galax- 
ies in GOODS-N by paper 1 in the table. Distributions of 
the best-fitting parameters are also shown in right panel of 
Figure [T3] The peak (i.e., mode) of the B/A distribution is 
determined from /i and a parameters, while /i^ and pa- 
rameters corresponds to a peak and width of the C /A dis- 
tribution, respectively. The single-zgso galaxies in GOODS-S 
show similar B /A distribution to those in the GOODS-N with 
a slightly larger width of GOODS-N sample (top right panel 
of Figure [T3]l. The best-fitting ji of the single-zgso galaxies 
in GOODS-S is -1.05^o 25, consistent with that in GOODS- 

N within Icr (/i = -0.95lo;i5)- Although a sHghtly larger 
best-fitting a in the GOODS-N samples, the a parameters of 
the samples in both fields are well in agreement with each 
other within 68% confidence level (Table|4|. The bottom right 
panel of Figure [13] shows the slightly different distributions 
of the intrinsic edge-on C/A ratio of the single-zgso galax- 
ies between in GOODS-S and in GOODS-N. However, their 
best-fitting and cr^, which correspond to C/A distribution, 
agree well to each other within Ict uncertainties. This indi- 
cates that the sBzK galaxies in both GOODS-S and GOODS- 
N are not round {B/A ~ 1.0) and seem to be flat (mean 
C /A = 0.25 — 0.28) when being observed in rest-frame UV 
wavelength. 

Since the intrinsic shapes of the single-zgso sBzK galaxies 
in both GOODS-S and GOODS-N fields are consistent with 
each other, we combined the samples together and re-derived 
the intrinsic shape in order to improve the statistics. The fit- 
ting results are shown in Table |4] The best-fitting parameters 
are similar to those determined in the separated fields but with 
smaller range of uncertainties. The results confirm that the 
sBzK galaxies that appear as a single component in the rest- 
frame UV wavelength show a flat (peak C/A — 0.26) and 
bar-like/oval shape (peak B/A = 0.71). 



6.3.2. Rest-Frame Optical in GOODS-S and SXDS 

Thanks to the publicly available data from CANDELS sur- 
vey, we are able to study the structure of sBzK galaxies at 
z 2 in the rest-frame optical wavelength. The apparent b/a 
histograms of the single-i/igo sBzK galaxies in GOODS-S 
and SXDS are shown in the left panel of Figure [14] The dis- 
tributions of samples in both fields look similar to each other; 
however, the K-S test indicates that these two histograms are 
drawn from different distribution. In any case, the distribu- 
tions roughly seem to be flat from b/a ^ 0.3 to b/a ~ 0.7 
and decline continuously toward b/a = 1.0. The distribu- 
tions of the single-iJi6o sBzK galaxies are still different from 
those of the local disk galaxies, which is also shown with the 
black solid line in the figure. Law et al. (2012) also studied 
the rest-frame optical morphologies of star-forming galaxies 
aX z ~ 1.5 — 3.6 and showed the b/a histogram that peaks at 
b/a 0.55. Their histogram seems to be similar to that of the 
sBzK sample in GOODS-S. 

We determined the intrinsic shape of the single-iJigo sBzK 
galaxies in the identical manner as done for the single- 
zg50 samples. The fitting results are summarized in Table H] 
The best-fitting parameters are also shown in the right panels 
of Figure [14] as B/A and C/A distributions. Despite statis- 
tically significant difference of the b/a distribution, the best- 
fitting i^L of the single-i/igo sBzK galaxies in both fields are 
the same at ^ = —1.30. Meanwhile, slightly different a pa- 
rameters result in different B/A distribution as seen in the top 
right panel of Figure [14] The best-fitting <t value of samples 
in GOODS-S is still consistent with that in SXDS within 68% 
uncertainty. The C/A distributions of the sBzK galaxies in 
both fields show almost the same peak with the same width, 
since they have comparable fi^ and cr^ values. 

Because the intrinsic shapes derived separately in GOODS- 
S and SXDS are well in agreement with each other within 
1(7 uncertainties, we combined the samples and re-derived the 
intrinsic shape to increase the statistical significance. Fitting 
results of the combined samples are listed in Table [4] The 
peak B/A and C/A ratios of the combined samples are 0.79 
and 0.27, respectively, indicating that the sBzK galaxies are 
indeed not round even when observing in the rest-frame opti- 
cal wavelength. The best-fitting parameters of the combined 
samples are identical to those of the single-ifieo galaxies in 
SXDS alone. This is presumably due to the fact that number 
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Table 4 

Best-Fitting Results of Intrinsic Shape for sBzK Galaxies 
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Figure 14. Left panel: Normalized distributions of the apparent axial ratio (b/a) of the single-Higo sBzK galaxies at 2 ~ 2 in the GOODS-S (red dashed 
histogram) and in SXDS (magenta solid histogram) fields. The best-fitting models are shown with soHd lines with corresponding colors. Distribution of the local 
disk galaxies by Padilla & Strauss (2008) is also shown with a thin solid histogram. Right panel: Distributions of the best-fitting intrinsic axial ratios: B/A 
distribution on the top right panel and C/A distribution on the bottom right panel with colors identical to the left panel. Y-axis is in arbitrary unit. 



of the single-_ffieo sBzK galaxies in SXDS is approximately 
three times larger than that in GOODS-S. 

6.3.3. Comparison between Rest-Frame UV and Optical 
Wavelengths 

In this section, we discuss the structural difference between 
in the rest-frame UV and in the rest-frame optical wavelength 
in terms of the apparent axial ratios (5/a). The b/a compar- 
ison of sBzK galaxies that appear as a single component in 
both wavelengths in the GOODS-S is shown in Figure [15] 
The left panel shows that most objects show larger h/a ratios 
in the rest-frame optical wavelength (i/ieo images) than those 
in the rest-frame UV (z850 images), though some show the 
same or even lower h/a ratios in the rest-optical images. The 
comparison of b/a histograms is illustrated in the right panel 
of the figure. It is seen that the b/a histogram is slightly dif- 
ferent at 6/a > 0.5 in the way that larger fraction of the sBzK 
galaxies in the rest-frame optical wavelength shows b/a ratios 
more than 0.5. The median fe/a value in the rest-frame optical 
images (b/a = 0.56) is also larger than that in the rest-frame 
UV {b/a = 0.51). This difference may result in the difference 
in the derived intrinsic shape of the sBzK galaxies. 

We investigate the difference of the intrinsic shape of sBzK 
galaxies between in the rest-frame UV and in the rest-frame 
optical wavelengths by using the combined samples in both 
wavelength ranges, i.e., the single-zsso galaxies in GOODS-S 



and GOODS-N for rest-frame UV and the single-iJieo galax- 
ies in GOODS-S and SXDS for rest-frame optical, to improve 
statistical significance. Their fe/a histograms are shown in the 
left panel o f FigurefTSlalong wi th those of the local disk galax- 
ies by Pad illa & StraussI (120081) . Similar to what we found in 
the GOODS-S, the b/a distribution of the sBzK galaxies ob- 
served in the rest-frame optical wavelength shows a slightly 
larger fraction at b/a > 0.5 than those in the rest-frame UV 
wavelength. Comparing the histograms with those of the lo- 
cal disks, we even see more difference of the histograms at 
b/a > 0.5. In paper I, they discussed that difference in b/a 
histogram results in difference in the intrinsic shape (Figure 
12 in paper I). The b/a histogram at 6/a > 0.5 affects the 
intrinsic face-on shape, mainly the /i parameters. Higher 
fraction at 6/a > 0.5 histogram indicates rounder face-on 
shape. The effect can be seen in the right panel of Figure 
[T6l the top right panel shows the intrinsic B/A distributions 
derived for the sBzK galaxies in the rest-frame UV and op- 
tical wavelengths. The B/A ratio of the single-zgso sBzK 
galaxies peaks at B/A = 0.65, which is different from that 
of the single-i7i6o sBzK galaxies more than la uncertainty. 
The same trend of dramatically increase at b/a = 0.1 — 0.3 
of the b/a histograms results in the similar i^i^ for the samples 
in both wavelength ranges, whereas some excess fraction at 
b/a = 0.3—0.4 in the rest-frame UV sample causes difference 
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Figure 15. Left panel: Comparison of appai'eiit axial ratios (b/a) of the sBzK galaxies in tfie GOODS-S tliat appear as a single component in both zgso and 
-ff l5o images. Right panel: Normalized b/a histograms of the single-^gso (blue dashed histogram) and the single-Higo sBzK galaxies (red solid histogram) in 
the GOODS-S. 



UV (GOODS-N+GOODS-S) 
Optical (GOODS-S+SXpS) 




Figure 16. Left panel: Normalized distribution of the observed axial ratio (b/a) of the combined sample of the single-^gso (GOODS-N+GOODS-S) and single- 
H \&0 (GOODS-S-l-SXDS) sBzK galaxies. The best-fitting models are shown with sohd lines with corresponding colors. Distribution of the local disk galaxies 
bv lPadilla & Strau"s3 120081) is also shown with a thin sohd histogram. Right panel: Distributions of the best-fitting intrinsic axial ratio: B/A distribution on the 
top right panel and C/A distribution on the bottom right panel with colors identical to the left panel. Y-axis is in arbitrary unit. 



in the derived cr-y parameters. Their best-fitting /x^ are similar 
at = 0.26 - 0.27 or peak C/A of 0.26 - 0.27, whereas 
wider C/A distribution can be seen in the single- TJigo galax- 
ies (bottom right panel) due to slightly larger . However, 
they are still in agreement within la confidence level. In con- 
clusion, as seen in the right panel of Figure[T6]and TablelH the 
sBzK galaxies appearing as a single component in rest-frame 
UV show just slightly different B /A or intrinsic face-on shape 
as compared to those in rest-frame optical wavelength with 
similar C /A or intrinsic edge-on distribution. 

6.4. Apparent h/ a Distribution and Intrinsic Shape of 
Subsamples 

6.4.1. Rest-Frame UV in GOODS-S and GOODS-N 

In this section, we discuss the dependence of apparent b/a 
distribution (if any) on vari ous pro perties of the sBzK galax- 
ies. As described in section l6.3.1l the b/a distribution of the 
single-Z85o sample in both fields are drawn from the same 
distribution and their best-fitting shape parameters agree well 
within 68% confidence level. Besides, the 2^50 images in both 
fields have almost the same depth. Therefore, we combined 
together the sample from both fields to increase the statis- 
tics. Figure [TT] shows b/a histograms of subsamples of the 
single-Z850 sBzK galaxies in GOODS-S and GOODS-N. As 



described in section 15.11 and also in paper 1, the Sersic index 
can be derived accurately down to zg.go = 24.5 mag in the 
Z850 images, but we used the single-zgso galaxies of which 
Z850 magnitudes reach zg^o = 26.0 mag in order to increase 
the number of samples. Thus we made a subsample of galax- 
ies with Z850 < 24.5 mag to investigate the shape difference 
as compare to the whole sample. Figure [TtI' a) shows the ap- 
parent b/a distributions of both zgso < 24.5 mag and the 
whole samples of single-zgso galaxies. The distributions of 
both samples are close to each other in the way that the his- 
togram of Z850 < 24.5 sample shows a peak at b/a 0.3 
and gradually declines until b/a = 1.0. We also derived the 
intrinsic shape parameters for subsamples and summarized in 
Table [5] The derived intrinsic shape for the brighter subsam- 
ple agrees well with that of the whole sample within 68% 
confidence level. Figure [17] (b) shows the very similar b/a 
histograms between the whole sample and 0.5 < n < 1.5 
subsample, while the histogram of the 1.5 <n< 2.5 subsam- 
ple is little fuzzy. Their best-fitting parameters in Tableware 
still in agreement within la uncertainty, suggesting that div- 
ing the single-Z850 galaxies into subsamples of 0.5 < n < 1.5 
and 1.5 < n < 2.5 does not affect the statistical results of 
the intrinsic shape. It is interesting to note here that sBzK 
subsample with n > 2.5 shows different b/a distribution. 



Intrinsic Shape of sBzK in GOODS-S and SXDS 



17 



Single-Zg5Q samples 




(d) 



(b) 



Single-Zg5Q samples 

0.5 < n< 1.5 

1.5<n<2.5 

2.5 <n 



0.2 0.4 0.6 0.8 



Single-z^ samples 
B-K<2.0 
B-K>2.0 



Apparent b/a 
1 0.15 I 



(e) 



0.2 0.4 0.6 O.E 

Apparent b/a 



(c) 



Slngle-ZDCQ samples 
9.0 <log(H3i3r)< 10.0 




a 



0.2 0.4 0.6 

Apparent b/a 



Single-Zg5Q samples ■ 
1.4<2< 1.9 ■ 
1.9<z<2.5 ■ 



0.2 0.4 0.6 

Apparent b/a 



Figure 17. Apparent b/a histograms of the combined sample of single-zgso sBzK subsamples in GOODS-S and GOODS-N. Solid histograms represent all 
single-2g5o sBzK galaxies with 0.5 < n < 2.5. 



The K-S test of two distributions between the whole sam- 
ple (0.5 < n < 2.5) and the n > 2.5 subsample results in 
P = 0.02, meaning that we can reject the null hypothesis that 
they are drawn from the same populations with 2% signifi- 
cance level. Their b/a histogram peaks at b/a ~ 0.6 and has 
a lower fraction at b/a < 0.5, resulting in larger disk thick- 
ness C/A but comparable face-on shape {B/A) as compared 
with n < 2.5 sample. 

Figure[T2](c) shows the comparison of b/a histograms for 
the lower-mass (9.0 < \og{Mstar/MQ) < 10.0) and higher- 
mass {log{Mstar / Mq) > 10.0) subsamples. The lower-mass 
subsample shows similar b/a distribution as compared to the 
whole sample; the derived parameters of intrinsic shape also 
agree well with each other. The histogram of the higher-mass 
subsample is slightly different from that of the whole sam- 
ple with P = 0.21 from the K-S test; however, their intrin- 
sic shape parameters are still in agreement within the 68% 
uncertainty. We also divided the single-zgso sBzK galaxies 
into subsamples according to their B — K colors. Histogram 
of B — K < 2.0 subsample is well in agreement with that 
of the whole sample as seen in Figure [TT] (d). Their de- 
rived intrinsic shapes are also consistent with each other. For 
the B — K > 2.0 subsample, the b/a histogram peaks at 
b/a ~ 0.40 — 0.45, which is larger than the b/a peak of the 
whole sample. This suggests a slightly higher C/A thick- 
ness as also seen in Table |5] though an error is large. Figure 
[TtI (e) shows the b/a histograms of the 1.4 < z < 1.9 and 
1.9 < z < 2.5 subsamples. The redshifts are binned so that 
both subsamples have comparable number of galaxies. His- 
tograms of both subsamples seem to look similar to those of 



the whole sample; their best-fitting parameters are also con- 
sistent within Icr uncertainty. 

6.4.2. Rest-Frame Optical in GOODS-S and SXDS 

Similar to the rest-frame UV wavelength, we combine the 
single-Fi6o sBzK galaxies from GOODS-S and SXDS to- 
gether in order to increase the statistical significance. Al- 
though the depth of the Hiqq image in GOODS-S is currently 
shallower than that in SXDS, the resulting intrinsic shapes 
derived separately i n both fields agree well with each other 
as shown in section 16.3.21 Thus combining the sample is un- 
likely to cause any false result or discussion. 

The single-i^ieo sBzK galaxies were divided into subsam- 
ples according to their various properties and their b/a his- 
tograms are shown in Figure [18] Figure [18] (a) shows the 
b/a histogram of the sBzK subsamples with Hiqq < 24.0 
mag, which is the magnitude limit for determining the ac- 
curate Sersic indices (section [STb . Their b/a histogram is 
remarkably identical to that of the whole sample and their de- 
rived intrinsic shape parameters are also the same (Table [5]), 
indicating that using the whole sample to increase the statistic 
accuracy does not affect the results on intrinsic shape. Figure 
[T8l(b) shows that the b/a histograms of the 0.5 < n < 1.5 and 
1.5 < n < 2.5 subsample have no significant difference from 
those of the whole sBzK sample. The fitting results in Table|5] 
also show an agreement within la uncertainties. In contrast, 
the n > 2.5 subsample shows lower fraction at b/a < 0.5 
and larger fraction at b/a > 0.5, similar to that seen in the 
rest-frame UV wavelength. According to Table [5] the peak 
B/A (0.83) and peak C/A (0.34) of the subsample are both 
different from those of the whole sample; the n > 2.5 sub- 
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Table 5 

Best-Fitting Parameters of Intrinsic Sliape for sBzK Subsamples 



Rest-Frame Wavelength 
(Observed Field) 


Subsamples 


A* 


(J 


M7 






All single- Z850 
^850 < 24.5 mag 


1 nn+0-30 
-l-00_o.3o 


45+0-10 
0-50-0.30 


26+0-02 
0.30_o 04 


n 055+0-015 

u.uoo_o 010 
0.030_o 020 


UV (GOODS-S&N) 


0.5 <n< 1.5 
1.5 < n < 2.5 
2.5 < n 


n st;+o i5 
-o.»a_o 30 

-t -1 ^+0 30 


o.5olS;t^ 

0.3510-60 


^ „ .4-0 08 
>-'-^^_0.02 
09+0.66 
'-'-■'^-0.08 


n 050+0-020 

U.UOU_o 020 
„_„-t_o 040 

0.050to-o« 

110+0-240 
U.liU_o 040 




9.0 < log(M*/M0) < 10.0 
10.0 < log(M*/M0) 


-o.«a_o 00 

^•^°-0.30 


0.50t0-30 
0.35lJ-^ 


n 00+0.02 

26+0-** 
'J-^°-0.04 


030+0-020 
u.uou_o 000 
070+0-220 
U.U(U_o 020 




B -K <2.0 
B - K>2.0 


„i 00+015 
-1 00+0-00 

i-UU_o 30 


0.50lO;?o 


26+0-02 
Q o(,+0.04 


050+0-020 
U.UOU_o 040 
070+0-020 
U.Um_o 020 




1.4 < z < 1.9 
1.9 < z < 2.5 


-1 00+0-00 

oc:+0.75 
-0.85_o.3o 


0-50ii;?° 

o,-+0.60 


n 90+0. 06 

0-^°-0.04 

24+0-20 
'J-^^-0.04 


n 050+0-040 

U.UOU_o 040 

0.050tO-J20 




All single-Z/ieo 
tl < z4.u mag 


-1 3D+0-iO 
i-JU_o.oo 

1 qn+ooo 


50+0- -^0 

U.OU_o.20 

n qc+o.ao 

U.JO_o.oo 


27+0-01 
^•^'-0.02 
9C+0.04 
0-2t>„o.oo 


0-080t°-.^^^ 

n 070+0-020 
0.070_o.ooo 


Optical (GOODS-S+SXDS) 


0.5 < n< 1.5 
1.5 < n < 2.5 
2.5 < n 


1 itr + O.OO 
-i.iO_o.i5 

-1 .T^+0.30 
-1 60+1-20 

i-OU_o 45 


I'-^'J-o.oo 

„ „„-f0.45 

0-501^.32 
0.35lJ:35 


26+0-02 
^-^"-0.02 

n or.+0.08 

0-30Io.o4 

S4+0-64 
U.04_o 04 


n 070+0-000 

U.U/U_o 000 

r\ r.r>r.+0.040 

0.0901OO2O 

iso+0-100 

U.lciU_o 100 




9.0 < log(M*/M0) < 10.0 
10.0 < log(M*/M0) 


1 ipr + O.OO 
— i.lO_0 15 

_-! 4^5+0-15 
l-^°-0.15 


0.50l«;i^ 


26+0-02 
'^■^"-0.02 

n 90+0.04 

0-^8-0.00 


070+0-000 

U.Um_o 000 
090+0-020 

u.uyu_o 020 




_B - A" < 2.0 
B - X > 2.0 


1-J^'^-O.OO 

-1 so+0-15 


0.50l«;l^ 

0.50l°:?o 


26+0-02 
'J-^O-0.02 

24+0-02 
'J-^^-0.02 


070+0-020 
U.U(U_o 000 

090+0-000 
u.uyu_o 020 




1.4 < z < 1.9 
1.9 < z < 2.5 


-1 so+0-15 

_1 -1^+0.15 


0.651°;?° 
0.201?,;go 


-^0+0-02 

U.dU_o 02 

24+0-0-1 
1-^^-0.02 


n 090+0-000 
u.uau_o 020 

070+0-040 
U.U<U_o 020 



sample statistically shows rounder and thicker shape than the 
n < 2.5 sample. Although the n > 2.5 subsample does not 
show the similar shape to the local elliptical gal axies (peak 
B/A = 0.98 and peak C/A = 0.43: lPadilla & Sia uss 2008]), 
the trend of being thicker in shape than the n < 2.5 subsample 
is the same as that in local universe in the way that local el- 
liptical galaxies show almost similar face-on ratio but thicker 
edge-on ratio than the local disk galaxies. 

The b/a distributions of mass-divided subsamples are 
shown in Figure [18] (c). Similar to the single- 
Z850 galaxies, the lower-mass subsample with 9.0 < 
\og{Mstar/MQ) < 10.0 shows similar 6/a distribution to the 
whole sBzK sample, meanwhile the larger-mass subsample 
with log{Mstar/MQ) > 10.0 shows slightly different dis- 
tribution, fewer fraction at b/a < 0.5 and larger fraction at 
b/a > 0.5. The derived parameters show that the larger-mass 
sample is slightly rounder and thicker than the lower-mass 
one; however, they are still within la uncertainty. Figure [TSl 
(d) illustrates the b/a histograms of the B — K subsamples. 
The b/a histograms of both B — K subsamples are similar 
to each other and similar to the whole sample; their derived 
parameters are also consistent within Icr uncertainty. Like- 
wise, both histograms of the redshift subsamples are not sig- 
nificantly different from those of the whole sample as seen in 
Figure[T8](e). Their derived intrinsic shapes are also in agree- 
ment with each other. 



7. DISCUSSION 



7.1. Intrinsic Shape Comparison with Disk Galaxies in the 
Local Universe 

The apparent b/a hist ogram of the local disk galaxies by 
iPadilla & Strauss! d2008h is shown in comparison purpose in 
the left panel of Figure [I6] As described in the previous sec- 
tion, the b/a histogram of the local disks is significantly dif- 
ferent from the sBzK galaxies both in rest-frame UV and op- 
tical wavelengths. The peak value of the derived B/A distri- 
bution is 0.79 ± 0.03 for the sBzK galaxies (combined sample 
in the rest-frame optical), while it is 0.98 ± 0.01 for the local 
disk galaxies, indicating that the local disks are significantly 
rounder than the sBzK galaxies. The p^ parameters of the lo- 
cal disks differs from that of the sBzK galaxies in both wave- 
lengths, suggesting that the local disks are statistically flatter 
than the sBzK galaxies at z ^ 2 (right panels of FigurefTSTl. It 
is implied that the population of galaxies showing round and 
flat disk seen in the local universe is unlikely to be exist yet at 
z - 2. 

At this redshift, the outermost isophotes we can reach in 
this study are ^ 27 — 28 and ^ 27 mag arcsec^^ for /-t^sso 
and HHieo^ respectively. Due to the cosmological dimming 
effect, Pz850 of ^ 27 — 28 mag arcsec^^ at 2 2 corre- 
sponds to /is ^ 23 — 24 mag arcsec"^ at z ~ after simple 
K correction by median Zgso — J color of the sample, while 
A^ffieo ^ 27 mag arcsec"^ corresponds to fiy ^ 23.5 mag 
arcsec"^ at z ~ 0. In contrast, the intrinsic sh ape study for 
local disk galaxies by IPadilla & StraussI ( 120081) could reach 
surface brightness of 25 mag arcsec"^. Although we could 
not reach the same surface brightness as observed in the lo- 
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(b) 



Single-H^gQ samples 
0.5 < n < 1.5 
1.5 < n<2.5 
2.5 <n 




n. 
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Apparent b/a 



Single-H.,QQ samples 
B-K<2.0 
B-K>2.0 
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Figure 18. Same as Figure [TTl but for the combined sample of single-Higo sBzK galaxies in GOODS-S and SXDS. 



(c) 



Single-HHRQ samples 
9.0 <log(M3i3r)< 10.0 
log(Msar)> 10-0 



0.4 0.6 
Apparent b/a 



Single-H.,gQ samples ■ 
1.4<z< 1.9 ■ 
1.9<z<2.5 





cal universe, it is possible that the surface brightness we see 
corresponds to ^ 25 mag arcsec"^ by conside ring ^ 1 mag 
lumin osity decrease from z ~ 1 to ^ '--^ (e.g. JScarlata et all 
l2007h . This surface brightness correspon ds to far beyond 
bar ends seen in the local universe (e .g., lOhta et all 119901 : 
iKuchinski e'tani2000l: lOhta et al.ll2007h . In fact, we trace the 
isophotes out to ~ 1" radius at z ^ 2, which corresponds to 
^ 8 kpc at z ~ 0. If the structural evolution was not drastic, 
we already trace the disk structure far enough. 

7.2. Why Bar-Like?: Possible Origins of the Bar-Like/oval 

Structure 

As discussed in the previous section, the sBzK galaxies at 
z 2 show the similar light profile, effective radius, and 
stellar surface mass density to those of the disk galaxies in the 
local universe, suggesting they are likely to be disk galaxies. 
However, we found that the sBzK galaxies (at least those with 
a single component) are statistically in bar-like or oval shape 
instead of a round disk seen in the local disk galaxies. This 
result rises an interesting question of why or how the sBzK 
galaxies become bar-like. It is worth noting here that the bar- 
like structure seen in sBzK galaxies is unlikely to be direct 
progenitors of bars in local barred galaxies, because a fraction 
of barred galaxies decreases with increasing r edshifts and it 
is only 10-20% at z ~ 0.8 (ISheth et al.ll2008h . Considering 
also the larger size as compared with a typical bar length in 
local universe, we regard the bar-like shape seen in the sBzK 
galaxies as a global structure of galaxies. 

As discussed in paper I, one possibility is bar instability 
when the fraction of disk mass a painst halo mass w ithin a 
disk radius exceeds the threshold (Ostriker & Peebles '1973b. 
Using A^-body simulations. .Athanassoula & Misiriotia (.20021) 



showed that different central concentration results in differ- 
ent morphology of the bars. Their massive disk model (MD 
model), where disk dominates the inner part of the galaxy, 
shows an agreement of the bar structure with the intrinsic 
shape of the sBzK galaxies. 

The ba r structure can also be rnade by galaxy interac- 
tion (e.g.,'Noguchi"1987;'Geri n et al.|[T990HMiwa & Noguchil 
il998; Berentzen et al. 2004). The galaxy is able to develop 
into a bar structure when perturbed by the tidal force of an- 
other galaxy. The galaxy loses its angular momentum due 
to tidal perturbation, resulting in elongating bar structure 
jBerentzen et al.l 120041) . Although this scenario seems to be 
less likely as the sBzK sample that we study is the single com- 
ponent, which shows no clear sign for strong close galaxy in- 
teraction, the bar-like structure generated thr ough the galaxy 
interaction can last f or at least a few Gyr ( iNoguchil 119961 : 
iBerentzen et al.l2004l) . Thus it is still possible that the bar-like 
structure we see in the sBzK galaxies are probably the results 
of the galaxy interaction in the past. The pass-by galaxy in- 
teraction (galaxy harassment) and continuous minor merging 
process are also likely to be responsible for the bar-like shape 
of the sBzK galaxies. Studying galaxy shapes in cosmological 
simulations is desirable to answer which of the above scenar- 
ios is the most likely process for causing the bar-like structure 
of star-forming galaxies at high redshift. 

7.3. Do We See A Thick Disk? 

The intrinsic thickness of C/A ~ 0.3 is close to a typ- 
ical value for thick disks in the present-day disk galaxies 
(e.g., lYoachim & Dale antoni 120061) . Age of the stellar pop- 
ulation composing the thick disk is old; in the Milky Way 
galaxy, the age of the stellar population in the thick disk is 
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older th an ^ 8 Gyr and m ost of them are between 10 and 
13 Gyr jReddy et alj|2006l) . Hence, there is a possibility that 
the structure we see is the progenitor of thick disks, though 
the radial scale length may be slightly smaller than that of 
present-day thick disk. The difference of peak thickness in the 
rest-frame UV and optical is very small, suggesting almost no 
age and/or metallicity gradient in the vertical st ructure. The 
thick disk may be heated up from a thin disk (e.g.. lOuinn et alj 
11993b . In the redshift range we studied there seems to be a 
subtle evolution of thickness as seen in Table |5] I-L^ = 0.28 
in 1.4 < 2 < 1.9 while it is 0.24 in 1.9 < z < 2.5 in 
the rest-frame UV and 0.30 in 1.4 < z < 1.9 and 0.24 in 
1.9 < z < 2.5 in the rest-frame optical wavelength. However, 
the difference is probably not significant when considering the 
large uncertainty. 

One problem for th e thick disk interpretat i on ma y be the 
stellar mass of disks. lYoachim & DalcantonI (120061) showed 
A/thick/A/thin = 0.53(^^/100 km s-^)-^-^ for the local disk 
galaxies. Our sample galaxies show a typical stellar mass of 
lO^°A/0 and a scale length of 2 kpc, giving a circular velocity 
of ~ 150 km s^^ if the stellar mass dominates the total mass. 
Then the expected stellar mass ratio between thick and thin 
components at z ^ would be smaller than ^ 0.2, i.e., the 
stellar mass of the thin disk at z should be larger than that 
of the sBzK galaxies at z ~ 2 with at least factor of five. If this 
is the case, the data points for z ^ 2 in Figure [12] should shift 
toward right by ^ 0.7 dex or larger as the galaxies evolve, 
provided that the d isk scale lengths do not evolve so much 
(iTrujillo et al.ll2006l) . Such a large shift can not reproduce the 
distribution of local disks in the figure. Thus the structure we 
see may partly include the thick disk, but would not totally 
be the thick disk. It should be worth noting that most of our 
sBzK galaxies are expected not to be evolving into present- 
day elliptical galaxies, because faint {K ^ 23 mag) sBzK 
galaxies are expected to reside in less massive dark halos and 
their counterparts at z ^ are exp ected to be disk gal axies 
according to a clustering amplitude dHavashi et al.ll20()7l) . 

8. SUMMARY 

We study the intrinsic structure of star-forming galaxies at 
z ~ 2 selected as sBzK based on their B — z and z~K colors. 
1028 and 29835 galaxies were selected down to Ks ~ 24.0 
mag in GOODS-South and SXDS fields, respectively. For 
these objects, we made SEDs and derived the photometric 
redshifts and stellar masses. In the GOODS-S, the structure 
of the sBzK galaxies was studied in both Z85o and i/ieo im- 
ages obtained with HST, which correspond to the rest-frame 
UV 3000A) and optical wavelengths 5300A), respec- 
tively, while only the rest-frame optical structure was studied 
in SXDS due to the availability of the HST images. The 
iJieo images cover only part of both GOODS-S and SXDS 
fields; therefore, only a fraction of sBzK galaxies (724 in 
GOODS-S and 2500 in SXDS) are used to study the rest - 
frame optical structure. Similar to paper I (I Yuma et al.ll201 ll) . 
majority of the sBzK galaxies appear as a single component in 
high-resolution images. 57% (583) of the sBzK galaxies were 
classified as a single component in zgso images in GOODS- 
S, while 86% (626) and 82% (2044) are classified as a sin- 
gle component in iJieo images in GOODS-S and SXDS, re- 
spectively. The larger fraction of single-component galaxies 
seen in -ffieo images is partially due to their poorer resolution; 
however, more than half of the single-iJigo galaxies seen as 
multiple components in the zgso images reveal the real mor- 
phological difference of galaxies between in the rest-frame 



UV and in the rest-frame optical wavelengths. 

Structural analysis was done for the single-component 
sBzK galaxies by fitting 2D light distributions of galaxies 
with a single Sersic profile by using GALFIT. Most of our 
samples show Sersic indices of n = 0.5 — 2.5, suggest- 
ing a disk-like structure. Regardless of the observed fields, 
the effective radii are in the ranges of 1.0 — 3.0 kpc and 
1.5 — 4.0 kpc in the zgso and i/igo images, respectively. The 
median effective radius is 1.89 kpc for the single-zg5o galax- 
ies in GOODS-S, while they are 2.52 kpc and 2.62 kpc for 
the single-i7ieo galaxies in GOODS-S and SXDS, respec- 
tively. Ratio between the effective radii in zgso and iJieo im- 
ages (re.opt/re.uv) is 1.3 3, consistent with that of BM/BX or 
LBGs at similar redshift jSwinbank et al.ll20I0l) . Size-stellar 
mass relation of the sBzK galaxies with 0.5 < n < 2.5 indi- 
cates that most sBzK galaxies have comparable stellar surface 
mass density to disk galaxies at z ^ — 1, again suggesting 
they are disk galaxies. 

Although the Sersic index, effective radius, and stellar sur- 
face mass density of the sBzK galaxies are comparable to 
those of the local disk galaxies, it is important to study their 
intrinsic structure to see whether disk galaxies with similar 
shape as observed today already exist at z 2. Distribution 
of apparent axial ratios (b/a) of the galaxies was used to de- 
rived the intrinsic shape of galaxies. By assuming a tri-axial 
model with axes A > B > C, we found that the sBzK galax- 
ies have peak face-on B/A ratio of 0.70 and peak edge-on 
C/A ratio of 0.25 - 0.28 in the rest-frame UV wavelength, 
where we also included the samples in the GOODS-N stud- 
ied in paper I to improve statistical accuracy. For the rest- 
frame optical wavelength, the peak B/A and C/A ratios are 
0.77 — 0.79 and 0.26 — 0.27, respectively. The results indicate 
that the sBzK galaxies have almost the same intrinsic shape 
in both rest-frame UV and optical wavelengths, though with 
slightly rounder in the rest-frame optical images. 

In order to examine the dependence of the intrinsic shape, 
we divided the sBzK galaxies into subsamples according to 
their properties. First, we determined the intrinsic shape of 
the brighter sample, whose structural parameters can be de- 
termined more reliably, and found that their derived intrinsic 
shape is consistent with the whole sample regardless of the 
rest-frame wavelength ranges. Likewise, dividing the sample 
into 0.5 < n < 1.5 and 1.5 < n < 2.5 does not affect the de- 
rived intrinsic shape. We, furthermore, divided the sBzK sam- 
ple into subsamples by considering their stellar mass, B — K 
color, and photometric redshifts. Except for the larger-mass 
subsample (log(A/sfar A/q) > 10.0) that shows very slightly 
rounder and thicker shape, the intrinsic shapes of other sub- 
samples are in agreement with those of the whole sBzK sam- 
ple. 

Comparing the intrinsic shape of the sBzK galaxies at z 2 
to those of loc al disk galaxies (peak i3M = 0.95, peak 
C/A = 0.21; lPadilla&Straussll2008h . we found that the 
sBzK galaxies are rather in bar-like or oval shape and slightly 
thicker than round and flat disks seen in the local universe. 
There are several possibilities explain the origin of the bar- 
like structure of the sBzK galaxies. First is the bar instabil- 
ity in a massive disk galaxy where disk mass dominate the 
central part of the galaxy. Interaction between galaxies is an- 
other possible process in which the galaxy can evolve into bar 
structure after influenced by the tidal perturbation from the 
pass-by galaxy. Even though no sign of interaction is seen in 
our single-component sample, we cannot reject the possibility 
of this galaxy interaction scenario as the resultant bar struc- 
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ture from the galaxy interaction can remain for a few more 
Gyr. Moreover, the continuous minor merge could also be 
the case. Study of structure evolution in cosmological simu- 
lations is desirable. As the intrinsic C jA ratio of the sBzK 
galaxies are comparable to the typical value for thick disks in 
the local disk galaxies, we may see progenitors of thick disks. 
However, if all the stellar mass of the sBzK galaxies at z ~ 2 
is that of thick disk, the expected stellar mass of total disk at 
z ~ would be very large, suggesting they are not the di- 
rect progenitor of the present-day thick disks, though a part of 
them may make the thick disks. 
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